Loop Dynamics in DNA Denaturation 
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The dynamics of a loop in DNA molecules at the denaturation transition is studied by scaling 
arguments and numerical simulations. The autocorrelation function of the state of complementary 
bases (either closed or open) is calculated. The long-time decay of the autocorrelation function is 
expressed in terms of the loop exponent c both for homopolymers and heteropolymers. This suggests 
an experimental method for measuring the exponent c using florescence correlation spectroscopy. 
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The thermodynamic properties of DNA near the ther- 
mal denaturation transition have been extensively stud- 
ied during the last few decades 0, • At low temper- 
atures a small fraction of the base pairs are unbound, 
forming loops of fluctuating lengths. These loops increase 
in size as the temperature is raised, until the denatura- 
tion transition is reached and the two strands separate. 
Experiments using uv absorption and specific heat mea- 
surements have yielded valuable information on equilib- 
rium properties of DNA 2}. Recently, single molecules 
techniques, most notably Fluorescence Correlation Spec- 
troscopy (FCS) have been used to study dynamical prop- 
erties such as the temporal behavior of loops |3| . 

The main theoretical approach for studying DNA de- 
naturation has been introduced by Poland and Scheraga 
(PS) |4( and was used to analyze the case of homopoly- 
mers. It was found that the dependence of the entropy 
of a loop on its length plays a dominant role in deter- 
mining the thermodynamic behavior near the transition. 
On general grounds one can argue that the entropy of 
a loop of length n takes the form S = fcs log(fi(n)), 
where Sl(n) ~ s n /n c is the number of loop configura- 
tions. Here s is a model-dependent constant and c is a 
universal exponent whose numerical value has been de- 
bated over the years and was found to depend on the 
degree in which excluded volume interactions are taken 
into account 0, 0, @- When excluded volume interac- 
tions both within a loop and between the loop and the 
rest of the chain are taken into account one finds c ~ 2.12 
0,@. This result, which predicts a first order denatura- 
tion transition, has been verified numerically 0. While 
numerical studies of the model with excluded volume in- 
teraction yield a clear first order transition |8|], a direct 
experimental measurement of c is rather difficult and has 
not been carried out so far. Theoretical studies of the 
case of a heteropolymers suggest that disorder makes the 
transition of order higher than two 0, . 

In this paper we analyze the loop dynamics at the de- 
naturation transition. The analysis suggests a method 
for measuring the exponent c. We focus on predictions 
for FCS studies Q. In these experiments one monitors 
the state of a base pair (whether it is open or closed) as 
a function of time. The measured quantity is the base 
pair autocorrelation function Ci(t) — (u,-(0)uj(t)) where 



ui (t) = 1 , is a variable which indicates if base pair i is 
open (1) or closed (0) at time t. By analyzing the loop dy- 
namics using a scaling approach and by direct modelling 
we express the temporal behavior of the autocorrelation 
function at the transition temperature in terms of the 
exponent c. 

Previous analyses of the loop dynamics have concen- 
trated mainly on the off critical region [TTl ll2T ] . In these 
analyses the loop is assumed to be in thermal equilibrium 
throughout its evolution. As discussed below in detail, 
the validity of this assumption is not obvious. In this 
work we test this crucial assumption and demonstrate 
that it is valid at the denaturation transition. 

To proceed we consider the dynamics of a single iso- 
lated loop. In this approach one ignores processes like 
merging of loops and the splitting of a large loop into 
two or more smaller ones. This may be justified by the 
fact that the cooperativity parameter, which controls the 
statistical weight of opening a new loop, is estimated to 
be rather small, <tq w 10 -4 0. Thus splitting a loop 
into two is unfavorable. Also, the average distance be- 
tween loops, which within the PS model is proportional 
to l/co, is large, making the independent loop approxi- 
mation plausible. A loop may change its size by closing 
or opening of base pairs at its two ends. It survives as 
long as its two ends do not meet. Let G(n, t) be the 
survival probability of a loop of initial length n for time 
t. The equilibrium autocorrelation function, measured in 
FCS experiments, is given by 



C(t) 



=i P eq {n)nG(n,t) 



(1) 



where for simplicity of notation we have dropped the site 
index i. Here P eq (n) is the probability of having a loop of 
length n in equilibrium. The additional n factor accounts 
for the fact that site i may be in any of the n sites of the 
initial loop. Note that we assume that site i remains 
open as long as the loop survives. This assumption does 
not affect the behavior of the autocorrelation function in 
the scaling limit. An interesting configuration is created 
when one end of the loop is forced to be on a particular 
site. In this case no n factor is needed in Eq. ^and the 
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autocorrelation function for the end of the loop reads 



c E (t) 



P eq (n)G(n,t). 



(2) 



Experimentally, this autocorrelation function may be 
measured by studying a molecule which is clamped at 
one end with a mismatch near this end. The autocorre- 
lation function near the mismatch site yields C E (t). 

In the following we analyze the cases of homogeneous 
and heterogenous DNA. We show that in the homo- 
geneous case the autocorrelation decays at large t as 
C(t) ~ t 1 -"/ 2 for c> 2 while it remains finite, C(t) = 1, 
for c < 2. On the other hand we find C E {t) ~ i( 1 - c )/ 2 
for c > 1. Our analysis of heteropolymers suggests that 
the disorder average of the autocorrelation function be- 
haves as C E (t) ~ (lni) 2 ~ 2c for 1 < c < 3/2 and as 
Ce ~ (logi) -1 for c > 3/2. Here the overline denotes an 
average over disorder. 

Consider first the case of a homopolymer. In this 



case it has been shown that P eq {n) 



The 



correlation length £ diverges at the transition yielding 
P eq (n) ~ nT c . In order to estimate the survival proba- 
bility of a loop of length n we consider the dynamics of a 
loop under the assumptions discussed above, where loops 
are non-interacting and they do not split into a number 
of smaller loops. Similar to 0, ^| we further assume 
that the loop is in a local thermal equilibrium at any 
given time during its evolution. The validity of this as- 
sumption will be discussed in detail below. The loop free 
energy is thus given by / oc n/£ + clnn where n is the 
loop size. Within the framework of the Fokker- Planck 
equation, the probability distribution of finding a loop of 
size n at time t, P(n, 1), is given by 



dP(n,t) _ D d 
dt dn 



c 
n 



d_ 

dn 



P(n,t) 



(3) 



where D is the diffusion constant in base pair units. 
Here we have taken the continuum limit and assumed 
the dynamics to be over-damped. This equation has to 
be solved with the boundary condition P(0, t) = and 
initial condition P(n, 0) = 5(n — no). The survival prob- 
ability of the loop is given by G(n ,t) = f Q dnP(n,t). 
Using standard techniques [R|] it can be shown that at 
the transition temperature = 0) the survival prob- 
ability obeys the scaling form G(no,t) — g(Dt/nf ) ) with 
z = 2. The asymptotic behavior of the scaling function 
for small and large values of the argument is 

g(x) ~ 1 for x <C 1 ; g(x) ~ x~~^~ for x 3> 1 . (4) 

Using these properties it is easy to calculate the long-time 
behavior of the autocorrelation function (Eq. 



C(t) 



1 for c < 2 

t l-c/2 for c > 2 _ 



(5) 



Thus, the asymptotic behavior of C{t) could in principle 
be used to measure the exponent c. In particular it can 



be used to distinguish between the case of a continuous 
transition (c < 2), where C(t) = 1, and a first order phase 
transition (c > 2), where C(t) decays to zero. Similar 
analysis for the edge autocorrelation function leads to 
C E (t) ~ t(!- c )/ 2 for c> 1. 

A central assumption in the above analysis is that the 
loop is at local equilibrium at any given time. A priori 
this is not necessarily a valid assumption. The typical 
time for the survival of a loop of length n scales as n 2 . 
On the other hand the relaxation time of a loop configu- 
ration is also expected to scale as n 2 when hydrodynamic 
interactions are ignored (to be discussed below) . Thus it 
is not clear that during the evolution of the loop it is 
in local equilibrium. Away from the transition point the 
loop size changes linearly in time and therefore the as- 
sumption of local equilibrium is clearly not valid. In the 
following we introduce a simple model for studying the 
loop dynamics where hydrodynamic interactions are ig- 
nored. We find strong evidence that the local equilibrium 
assumption holds asymptotically even in this case. 




FIG. 1: A typical microscopic configuration of the loop in 
the RSOS model. Dashed lines indicate possible dynamical 
moves of the interface. 

To this end we introduce and analyze a simple model 
for the loop dynamics. This dynamics is described by a 
fluctuating interface interacting with an attractive sub- 
strate in d = 1 + 1 dimensions. Here the interface height 
variable corresponds to the distance between complemen- 
tary bases. The interface configurations are those of a 
restricted solid on solid (RSOS) model defined as follows 
(see Fig. 0: Let hi = 0, 1, 2 . . . be the height of the in- 
terface at site i. The heights satisfy \hi — h i+ i \ = ±1 and 
hi > 0. Consider a loop between sites and n (where n 
is even) as shown in Fig. ^ Outside the loop the inter- 
face is bound to the substrate so that h-2k — h n+ 2k = 
for k = 0, 1, . . . while for < k < n the height hi can 
take any value which is consistent with the RSOS con- 
ditions and is non-negative. We consider a dynamics in 
which the loop is free to fluctuate and one of its ends 
(n > 0) is free to move with the following rules. For sites 
2 < i < n — 1 the height is updates according to 



hi ± 2 



(6) 



with rate 1 as long as hi > and the RSOS condition is 
satisfied. For i — n the loop length is changed according 
to the following rules 



n — > n + 2 with rate a/4 
n — > n — 2 with rate a , 



(7) 
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where n can decrease only if h n -2 = 0. 

In principle one should let both ends fluctuate. How- 
ever, for simplicity, we consider the case where one of the 
ends is fixed. It is straightforward to verify that the num- 
ber of configurations of a loop of size n is given by 2 n /n c 
with c = 3/2 for large n. Thus this model corresponds to 
a particular value of c. However we expect similar results 
to hold for other values of c as well. The ratio between 
the two length changing processes is chosen such that in 
the large n limit the loop is not biased to either increase 
or decrease. This corresponds to the model being at the 
denaturation transition point. The parameter a deter- 
mines the rate of the length changing processes: a — 
corresponds to the dynamics of a loop of fixed length and 
as a is increased the length changing processes become 
faster. In a realization of this dynamics at any given step 
one of n possible moves is chosen. Of these, n — 2 moves 
correspond to an attempted update of the height at sites 
2, 3, . . . , n — 1. The other two moves correspond to an 
attempt to update the edge by a move either to the right 
or to the left. One attempted move of the edge defines 
a Monte Carlo sweep. On average this amounts to up- 
dating all sites every two sweeps. The numerical studies 
described below are done using a = 1. 

In order to test the validity of the Fokker-Planck equa- 
tion J2J for describing the dynamics of a loop we simu- 
lated the dynamics of the model and calculated the sur- 
vival probability of a loop of initial length no- To this 
end an initial configuration of a loop of fixed length hq is 
generated with the correct equilibrium weight. Starting 
from this initial configuration the dynamics is carried out. 
The results are summarized in Fig. |2 where the survival 
probability is plotted as a function of the scaling variable 
t/n,Q for several values of the loop size ran. A very good 
agreement with the predicted survival time obtained from 
the solution of the Fokker-Planck equation © is found. 
However the optimal data collapse takes place atz« 2.2 
rather than z = 2. If this value of z remains valid in the 
limit of large n it would imply that the local equilibrium 
assumption is not valid. 

In the following we argue that the value z = 2.2 is due 
to finite size effects and we expect that for large systems 
the value z = 2 is recovered. To this end we calculate 
numerically the variance of the loop size 



w 2 (t) = ((n(t) - (n(t))) 2 ) , 



(8) 



where (•} denotes an average over realizations of the dy- 
namics. We show that asymptotically it grows linearly 
with time. This result indicates that the dynamical ex- 
ponent is in fact z = 2 and that the deviations we observe 
for small uq are due to finite size effects. We proceed by 
defining a variable <J+(t) which takes the value 1 if the 
length of the loop increases at time t and otherwise. 
Similarly, we define <J-(t) and ao(t) for steps which de- 
crease the loop size and steps in which the loop size does 
not change, respectively. Clearly o+(t) + <r_ (t) + (Jq (t) = 
1. The dynamics of the chain Q, implies that in the 
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FIG. 2: Data collapse of the survival probability (averaged 
over 4 ■ 10 4 realizations) for some values of no with dynamical 
exponent z — 2.2. The line corresponds a numerical solution 
of Eq. 



limit of large n one has 

(a+(t)) = (£T_(i)> - a/8 ; (a (t)) = 1 - a/4 , (9) 

where a = a/ max{l,a}. Denoting U(t) = a+(t)—a-(t), 
it is easy to see that 



Aw 2 {t) 



At 



w 2 (t) - w 2 (t - 1) 



= 4(1^) 2 }+8]T([/(t)[/(;)}, (10) 



T=l 



where 



(U(r)U(t)) = (a + (r)a + (t)) + (<7_ (?>-(*)> 

- (a_(r)<7 + (i))-(<7 + (r)<7_(t)). (11) 

It is evident that a loop increasing step at time t, 
(<r+(t) = 1), is uncorrelated with steps which took place 
at time r < t. Thus (o- + (-r)er + (i)) = (er_ (t)(7 + (t)} = 
a 2 /64. Numerically we find (cr_(r)cr_(f)) = a 2 /64 (see 
Fig. |2J). Using these result we finally obtain 



Aw 2 (t) 
At 



(12) 



with (er + (T)(7_(i)} c = (er + (T)<7_(t)} — Numerical sim- 
ulations of the dynamics show strong correlation between 
o+(t) and (T-(t) with an algebraic decay in t — r (see Fig. 
[3J|. It is interesting to note that the dynamics of the chain 
induces such long range temporal correlations mediated 
by the loop dynamics. 

By extrapolating the sum on the right hand side of 
Eq. E| using the asymptotic form At' 1 with A ks 0.015 
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the length changing rate ai) are taken to be uncorrelated 
quenched random variables. As the edge moves m steps 
the binding energy changes by A(m) = Y^iLi Vify- The 
rate oti is simply related to the binding e nerg y r](i), sim- 
ilar to the dynamics of DNA unzipping [14(. Since the 
variance of A(m) grows linearly with m we can safely 
neglect the effect of logarithmic correction to the loop 
entropy on the dynamics. At the transition point the 
dynamics of the loop length corresponds to that of an 
unbiased walker on a random forcing energy landscapes. 
It is known that the probability of a walker (representing 
the edge of the loop) not to return to the origin on such 
a disordered energy landscape, GAn.t), has the scaling 
form Gd(n, t) = g d ((logi) 2 /no) The asymptotic be- 
havior of the scaling function for small and large values 
of the argument is 



FIG. 3: Correlation functions of the a variables as obtained 
by averaging over 1.9 ■ 10 J realizations. 



and 7 ps 1.2, deduced from Fig. [21 we find that the sum 
converges to a value « 0.84 < a = 1 indicating that 
w 2 (t) 0.16i at large t, which in turn yields z = 2 
(see FigQJ. The slow power-law convergence towards 
the asymptotic value implies that it may require large 
systems to observe the long time behavior of the auto- 
correlation function, (JSJ). 



g d {x) - 1 for x < 1 ; g d (x) 



for x > 1 . (13) 
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FIG. 4: Temporal behavior of Aw 2 / At showing its slow con- 
vergence to its asymptotic non-vanishing value. 

We now turn to the heteropolymer case. To avoid 
the typical denominator problem in disordered systems 
we only consider the autocorrelation function (7^(4) and 
study its disorder average. For a heteropolymer the bind- 
ing energy of the zth base-pair, rj(i), (and therefore also 



This result is universal and independent of the potential 
realization. Note that in this case there is a separation 
of time scales where the typical survival time of the loop 
is much longer than the loop relaxation time. Thus the 
use of local equilibrium dynamics is clearly justified. 

To complete this analysis one has to evaluate the equi- 
librium loop size distribution. Extensive numerical stud- 
ies suggest that the disorder average loop statistics re- 
main of the same form with the same exponent c as in 
the case of homopolymers 0, ■ 

Combining this with the universal form of the survival 
probability we finally reach the asymptotic form of dis- 
ordered average autocorrelation function: 



C E (t) 



(logi) 2 - 2c 
(log*)" 1 



for 
for 



1 < c < 3/2 
c > 3/2 



(14) 



We conclude with a comment on hydrodynamic inter- 
actions. In the present study these interactions have not 
been included. It is well known that relaxation processes 
of polymers in solutions are faster when hydrodynamic 
interactions are taken into account. Within the Zimm 
model it scales with the length of the polymer as 
r? v . Recent experiments on single-stranded DNA have 
measures a scaling n ' This fast relaxation time 

lends further support for the local equilibrium hypothe- 
sis applied in this work. 
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